The Shear Mode of Multi-Layer Graphene 
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We uncover the interlayer shear mode of multi-layer graphene samples, ranging from bilayer- 
graphene (BLG) to bulk graphite, and show that the corresponding Raman peak measures the 
interlayer coupling. This peak scales from^43cm -1 in bulk graphite to^31cm -1 in BLG. Its low 
energy makes it a probe of near-Dirac point quasi-particles, with a Breit-Wigner-Fano lineshape due 
to resonance with electronic transitions. Similar shear modes are expected in all layered materials, 
providing a direct probe of interlayer interactions. 



Single Layer Graphene (SLG) has high mobility and 
optical transparency, in addition to flexibility, robust- 
ness and environmental stability [1, 2]. As the knowledge 
of the basic properties of SLG increases, an ever grow- 
ing effort is being devoted to a deeper understanding of 
Few Layer Graphene (FLG) samples [3-5], and to their 
application in useful devices. For example, since SLG 
absorbs 2.3% of the incident light [6], FLG can be used to 
beat the transmittance of Indium Tin Oxide(~90%)[2], 
and to engineer near-market transparent conductors [7], 
exploiting the lower sheet resistance afforded by combin- 
ing more than one SLG [2, 7]. Bilayer graphene (BLG) is 
a tunable band gap semiconductor [8] , tri-layer graphene 
(TLG) has a unique electronic structure consisting, in 
the simplest approximation, of massless SLG and mas- 
sive BLG subbands[9-ll]. FLG with less than 10 layers 
do each show a distinctive band structure [11]. The lay- 
ers can be stacked as in graphite, or have any orientation. 
This gives rise to a wealth of electronic properties, such 
as the appearance of a Dirac spectrum even in FLG [12]. 

There is thus an increasing interest in the physics and 
applications of FLG. Optical microscopy can count the 
number of layers[13, 14], but does not offer the insights 
of Raman spectroscopy, being this sensitive to quasi- 
particle interactions [15]. Raman spectroscopy is one of 
the most useful and versatile tools to probe graphene 
samples[15, 16]. The measurement of the SLG, BLG, 
and FLG Raman spectra [15] triggered a huge effort to 
understand phonons, electron-phonon, magneto-phonon 
and electron-electron interactions, and the influence on 
the Raman process of number and orientation of layers, 
electric or magnetic fields, strain, doping, disorder, edges, 
and functional groups [16]. 

The SLG phonon dispersions comprise three acous- 
tic and three optical branches. A necessary, but not 
sufficient, condition for a phonon mode to be Raman 
active is to satisfy the Raman fundamental selection 
rule, i.e. to be at the Brillouin Zone centre, T, with 
wavevector q « 0[17]. SLG has six normal modes at T: 
A2 U + B2 g + Ei u + E2 9 [18]. There are two degenerate 



in-plane optical modes, E^, and one out-of-plane opti- 
cal mode £>2#[18]. E2 g modes are Raman active, while 
B2 g is neither Raman nor IR active [18]. In the case of 
graphite there are 4 atoms per unit cell, and only half of 
them have fourth neighbors that either lie directly above 
or below in adjacent layers. Therefore the two atoms of 
the unit cell in each layer are now inequivalent. This dou- 
bles the number of optical modes, and is responsible for 
the IR activity of graphite[18]. All SLG optical modes 
become Davydov-doublets in graphite: E<i g generates a 
IR active E\ u and a Raman active E2 gi E>2 g goes into an 
IR- active ^2 n , and an inactive B2 9 . The zone boundary 
acoustic modes fold back to the zone centre as rigid layer 
modes: an optically inactive B2 g and a Raman active 
E2 g . The acoustic modes remain and Ei n [18]. Thus 
for graphite[18, 19] T = 2(A 2u + B 2g + E lu + E 2g ). There 
are now two Raman active E 2g modes, each doubly de- 
generate. The high frequency E2 g mode is responsible 
for the well-known G peak, measured and discussed in 
thousands of papers to date for any carbon allotrope[20]. 

Here we focus on the low energy E2 g mode. This is a 
doubly degenerate rigid layer shear mode, involving the 
relative motion of atoms in adjacent planes. It was first 
measured in 1975 by Nemanich et al.[21] in bulk graphite 
at~42cm _1 . We uncover the equivalent mode for FLGs, 
and show that it provides a direct measurement of the 
interlayer coupling. For this reason we name C the corre- 
sponding Raman peak. On one hand, the C peak energy, 
E(C)~5meV, is much lower than the notch and edge fil- 
ter cuts of most Raman spectrometers, and its intensity 
is much smaller than the G peak. This explains why it 
was not seen thus far in FLG and, even for graphite, it 
was reported only in an handful of papers [2 1-23], with 
no firm agreement on position and width. On the other 
hand, this makes it a probe of the quasi particles near 
the Dirac point. 

Raman measurements are performed at 633nm in 
backscattering geometry using a Jobin-Yvon HR800 sys- 
tem. The laser plasma lines are removed using a band- 
pass filter, since those would appear in the same spectral 
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FIG. 1. (al) Unpolarized and polarized Raman spectra of 
3LG on SiO2/Si(001) in the C peak region and (a2) in the 
G peak region. (b) Optical micrograph of FLG sample. 2LG- 
A/11LG-A and 2LG-B/11LG-B denote supported and sus- 
pended flakes, respectively. (cl) S/AS spectra of supported 
samples in the C peak region. (c2) S-spectra of supported sam- 
ples in the G/2D peaks region. (dl) S/AS spectra of suspended 
samples in the C peak region. (d2) S-spectra of suspended sam- 
ples in the G/2D peak region 



range as the C peak. A typical laser power of 0.5mW 
is used to avoid sample heating. Detection of Raman 
modes down to^lOcm -1 is possible by using BragGrate 
notch niters with optical density 3, and with a spectral 
bandwidth^5-10cm _1 . Three filters for each excitation 
are necessary to suppress the Rayleigh signal, which is 
typically 10 9 -10 12 stronger than the Raman signal[17]. 
Ar gas is flown on the sample to remove the low- frequency 
Raman modes from the air. We use a 100 x objective with 
NA=0.90. A 1800 lines/mm grating enables us to have 
each pixel of the charge-coupled detector cover~0.4cm _1 . 
A spectral resolution^. 5cm -1 is estimated from the 
width of the Rayleigh peak. 

The easiest way to get high quality SLG and FLG 
is by graphite exfoliation [24] on Si02/Si, to enhance 
visibility [13, 14]. Often the Si is doped, to be used as back 
gate [24]. However, this poses a problem for low frequency 
Raman measurements. The incident light can excite car- 
riers in doped Si, producing a strong background [25], 
that can overshadow the signal of FLG with less than 6 
layers. One approach to overcome this issue is to perform 
polarized Raman measurements, since this background is 
strongly suppressed in cross polarization [25]. Fig.l(al) 
shows the unpolarized Raman spectrum (top graph) of 



3LG on SiO2/Si(001), as well as the polarized one (bot- 
tom graph) with incident light along [1-10] and scattered 
light analyzed along [110]. A large substrate background 
is observed in the unpolarized measurement, while in the 
polarized one the Si mode and its low-frequency back- 
ground are suppressed, thus revealing a peak~37cm _1 . 
However, polarized Raman yields low overall intensity, 
and affects the 2D to G ratio. Thus, in order to detect 
the C peak in unpolarized measurements, we take a differ- 
ent avenue. We use low doping Si (resistivity>2000^.cm) 
and suspend the FLG on~2-5/mi etched holes, as shown 
in Fig. lb. The number of layers is identified by a com- 
bination of 2D peak Raman spectroscopy [15] and optical 
contrast on the supported section of the flake [13, 14]. 
Figs.l(cl,c2) plot the Stokes (S) and Anti-Stokes (AS) 
spectra for supported and suspended BLG, 11LG and 
bulk graphite. We use the notation NLG to indicate FLG 
with N layers. Thus 1LG=SLG; 2LG=BLG, 3LG=TLG, 
while, e.g., 11LG means 11 layers. In the suspended 
flakes the C peak is clearly seen. On the other hand, 
the supported ones show the Si background, Fig.l(cl). 
While for 11LG and bulk this does not overshadow the C 
peak, for fewer layers this covers the C peak, to the point 
that for supported BLG we cannot detect any C peak. 

Given the low E(C), the S/AS intensity ratio is close 
to 1. Similar to the G peak [26], but unlike D and 
2D [26], the S/AS C peaks are symmetric relative to the 
Rayleigh line. This allows us to precisely determine 
the C peak position, Pos(C), as [Pos(C)5-+Pos(C)as']/2. 
We get Pos(C)-31cm- 1 for BLG;-42.7cm- 1 for 11LG 
and~43.5cm _1 in bulk graphite. By assuming a 
lorentzian lineshape, we derive the Full Width at Half 
Maximum, FWHM, to be FWHM(C)-1.2cm- 1 . Consid- 
ering the^0.5cm _1 spectral broadening of our spectrom- 
eter, we derive an intrinsic linewidth~0.7cm _1 . From the 
S/AS ratio we estimate the local T on the sample as [17], 
T = hu/k B ln{I(C)s/I(C) AS • {[ujl + Pos(C)]/[uj l - 
Pos(C)]} 4 }, where h is the reduced Planck constant, 
ujl is the laser frequency, ks is Boltzmann constant, 
I(C) s /I(C)as is the C peak S/AS intensity ratio. This 
gives T~300K, indicating negligible laser heating. 

Figs.2(al,a2) plot the Raman spectra for a set of sam- 
ples of increasing thickness. Fig. 2(b) shows the fitted 
Pos(G) and Pos(C) as a function of 1/N, where N is 
the number of layers. While Pos(G)~1581cm _1 with no 
significant change with N, Pos(C) increases from BLG 
to bulk graphite. Note that in Fig. 2(b) the spectral 
range used to plot the G and C peak data is the same 
(16.5cm _1 ), thus the C peak shift with N is truly repre- 
sentative of its much stronger variation when compared 
to G. The ratio of peak heights at 633nm, I(C)/I(G) 
is-0.01, 0.06, 0.06 for BLG, 11LG and bulk graphite, 
while that of integrated peak areas, A(C)/A(G)~0.001, 
0.006, 0.007. These change with excitation energy, e.g., 
A(C)/A(G)~0.0014 for bulk graphite at 532nm. Since 
these ratios depend on the Electron Phonon Coupling 
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FIG. 2. (al) S/AS C peak spectral region. (a2) G peak spectral region. (b) (solid black circles) Pos(G), (open blue circles) 
Pos(C), as a function of inverse layer number. The red dash-dotted line is a plot of Eq.3, open diamonds are DFT calculations 



(EPC), this immediately indicates that EPC(C) is much 
smaller than EPC(G). 

The Pos(C) dependence on the number of layers can be 
explained considering a simple linear-chain model. For 
FLG with N layers, there are 2N atoms per unit cell. 
The corresponding in-plane optical modes consist of N 
degenerate pairs of in-plane stretching modes, and N-l 
degenerate pairs of in-plane shear modes between neigh- 
boring layers. We assume that a layer interacts strongly 
only with adjacent layers and that the strength of this in- 
terlayer coupling is characterized by an inter-layer force 
constant per unit area, a. The N-l shear modes of a 
NLG can be computed by diagonalizing the correspond- 
ing NxN (tridiagonal) dynamical matrix. The frequency 
uji (in Hz) of the i-th vibrational mode is given by: 



u I 



cos 
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where i=2,...N. /i=7.6x 10" 27 Kg/A 2 is the SLG mass per 
unit area. The corresponding z-th displacement eigenvec- 



tor is given by: 
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where j labels the layers. The highest frequency mode 
(for i=N) is Raman- active. Here adjacent layers move 
out-of-phase in the direction parallel to the planes. In the 
case of graphite N —> oo and = Pos (C)oo x (2ttc) = 



2y / a//i, where c is the speed of light in cm/s. This is the 
doubly degenerate E 2g shear mode responsible for the C 
peak, see Fig. 3(a). Thus, Pos(C)at (in cm -1 ) for a NLG 
is given by Eq.l setting i=N: 



(3) 



In BLG, N=2, and Pos(C) 2 = V^y^a/Jt x JL, i.e. V2 
smaller than Pos(C)oo, corresponding to bulk graphite, 
in excellent agreement with the experiments. In fact, the 
dash-dotted line in Fig. 2(b) shows that Eq.3 describes all 
the experimental data, thus validating our simple model. 
The only unknown parameter in Eq.3 is the inter layer 
coupling strength. By fitting the experimental data we 
can directly measure it. We get a ~ 12.8xl0 18 N/m 3 . 
Thus, in Bernal stacked FLG, the C mode hardening is 
not due to a variation of interlayer coupling, but rather to 
an increase of the overall restoring force going from BLG 
to bulk graphite. For a given N, we expect variations 
of Pos(C) if the interlayer coupling is modified, e.g. by 
changing the spacing or relative layer orientation (in the 
latter case we also expect mode splitting). 

These results are further confirmed by density func- 
tional theory (DFT) calculations. These are performed 
using DFT and density-functional perturbation (DFPT) 
theory as implemented in the PWSCF package of the 
QUANTUM-ESPRESSO distribution, within the Local 
Density Approximation(LDA), and ultrasoft pseudopo- 
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FIG. 3. (a) The two degenerate E23 shear modes in graphite. (b) Symmetry and ab-initio frequencies for each shear mode. The 
Raman and IR active modes are identified. 



tentials generated using the RRKJ approach. The cutoffs 
are 40Ry for the wave functions, and 480Ry for the charge 
density. The Brillouin zone is sampled on a 42 x 42 x 16 
Monkhorst-Pack mesh for bulk graphite and 42 x42 x 1 for 
SLG and BLG. NLG are modeled using supercell configu- 
rations, with periodic replicas separated by 10 A vacuum 
in the perpendicular direction. The electron-phonon and 
phonon-phonon matrix elements, as well as the anhar- 
monic contribution to the C mode linewidth, are com- 
puted using the approach of Ref.[27]. The EPC contri- 
bution to the linewidth is computed using an interpola- 
tion based on maximally-localized Wannier functions as 
implemented in the EPW code [28]. This is a computa- 
tionally efficient approach allowing very fine sampling of 
the Brillouin zone (meshes of several million points are 
needed to get accurate phonon linewidths). 

Fig. 3(b) plots the in-plane shear modes for 2-5LG and 
bulk graphite. For a given N, there are N-l shear modes, 
either Raman or IR active or both, but, for N>2, N-2 
of those have a different displacement pattern compared 
to C, since not all the neighboring layers vibrate out-of- 



phase. The highest frequency Raman active mode corre- 
sponds to the C peak. We expect the other Raman active 
modes to have a much weaker intensity compared to the 
C peak, as a result of a smaller EPC, also confirmed by 
DFT. More work is needed to detect those modes. 

We get an excellent agreement between our DFT fre- 
quencies and the experimental data, as indicated by 
the open symbols in Fig. 2(b). This might seem sur- 
prising, since local or semi-local exchange correlation 
functionals may not properly describe Van Der Waals 
interactions [29]. However, it was shown that in bulk 
graphite all phonon dispersions are well described by 
DFT, both in LDA and in the generalized gradient ap- 
proximation (GGA), even in absence of Van der Waals 
(VdW) interactions in the functional, provided that the 
correct geometry (i.e. interlayer spacing) is used [30]. 
This occurs because VdW interactions can give a signifi- 
cant contribution to the total energy (hence determining 
the c/a ratio, where c is the interlayer spacing and a the 
in plane lattice constant), but give a negligible contribu- 
tion to second derivatives (i.e. the phonons). Hence, if 



5 



the correct geometry is used, phonon dispersions are well 
reproduced. LDA does provide excellent geometries: c/a 
for graphite is 2.725, 2.74, 2.91, for experiments, LDA, 
VdW-DFT, respectively. For BLG, LDA and VdW-DFT 
give c/a=2.74, 2.90, i.e. VdW-DFT predicts the same 
spacing in bulk graphite and BLG, while LDA consis- 
tently predicts a smaller value than VdW-DFT, but in 
excellent agreement with experiments. This is confirmed, 
independently, for bulk graphite and any FLG, by the 
very good agreement between our LDA calculations and 
the measured FLG C modes. 

We now consider FWHM(C). Two factors contribute 
to the linewidth of the Fi2 g Raman modes in graphene 
and graphite: the EPC term[31, 32] and anharmonic 
phonon-phonon interactions [27]. In the absence of dop- 
ing, FWHM(G)-12-14cm- 1 in SLG and bulk graphite, 
mostly due to the dominant EPC contribution [31, 32], 
the phonon-phonon one being— 1.7cm _1 [27]. The experi- 
mental FWHM(C) is much smaller not just compared to 
the overall FWHM(G), but also with respect to the non- 
EPC component of FWHM(G). This immediately indi- 
cates a much smaller EPC(C) than EPC(G), consistent 
with the much smaller C peak intensity. Our DFT cal- 
culations give FWHM(C)-0.3cm- 1 at 300K in graphite, 
in reasonable agreement with experiments. The EPC 
contribution being— 0.05cm _1 , and the phonon-phonon 
one— 0.25cm _1 . The anharmonic term consists of three- 
phonon decay (30% of the total anharmonic linewidth at 
300K) and absorption (70%) processes. The shear mode 
splits mainly into two out-of-plane ZA bending modes, 
at q and -q, close to T. The absorption processes are 
dominated by the merging of the shear mode and a ZA 
mode into an out-of-plane ZO' bending mode (the prime 
indicates an optical mode where the two atoms in each 
layer of the unit cell of graphite vibrate together, but 
out-of-phase with respect to the two atoms of the other 
layer), see Fig. 4a. We expect the anharmonic linewidth 
not to change significantly in NLG, since the available 
phase space of decay /absorption channels in these sys- 
tems is very similar to graphite. Also, our calculations 
for BLG, TLG and 4LG show that the EPC contribu- 
tion to the linewidth is nearly independent of the num- 
ber of layers. Thus, DFT indicates the overall FWHM 
not to change significantly with N, in agreement with 
experiments. Note that, if we take a 4LG as an exam- 
ple, the EPC contribution to FWHM(G) is-150 times 
bigger than for FWHM(C). In turn, this is— 15 times 
bigger than the EPC contribution to the other Raman 
active shear mode— 17cm _1 , confirming the expectation 
that the other C modes would be challenging to detect. 

We now examine more closely the C peak shape. 
Fig. 4(b) shows that it can be well fitted with a Breit- 
Wagner-Fano (BWF). In general, this arises as quan- 
tum interference between a Raman allowed phonon and a 
continuum of Raman active electronic (or multiphonon) 
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FIG. 4. (a) Schematic representation of the anharmonic decay 
channels for the C mode in bulk graphite, (b) C peak for 3LG 
and bulk graphite, with BWF flt.(c) Schematic band structure 
of 3LG close to K. The gray regions highlight transitions near 
K that could resonate with the C mode. Red arrows indicate 
transitions with the same energy as the C mode, E(C). 



transitions [33]. The BWF lineshape is [33]: 

[i + 2(^-^o)/(gr)] 2 



[i + 4(^-^ ) 2 /r 2 ] 



(4) 



where Io, c^o, T and l/\q\ are the intensity, uncoupled 
mode frequency, broadening parameter and coupling co- 
efficient. The peak maximum is at (jO max =ujQ + T/2q, 
while its FWHM=r(<? 2 + l)/\q 2 - 1|. In the limit 
a Lorentzian lineshape is recovered, with FWHM=r and 
u m ax=uo. Pos(C) in the BWF fit is T/2\q\ (-0.3cm- 1 ) 
higher than in a Lorenzian fit. 

We find a smaller l/\q\ when we use a laser power high 
enough to shift the G peak, i.e. to heat the sample. 
Thus, in our low power experiments, the possible laser- 
induced electron-hole plasma is not the cause of the ob- 
served BWF lineshape. We also find that the C mode of 
bulk graphite at 77K has the same q as at room temper- 
ature, in contrast to what expected if the BWF would be 
due to a multiphonon resonance [34]. We thus attribute 
the BWF lineshape to quantum interference between the 
C mode and a continuum of electronic transitions near 
the K point. The band structure of Bernal-stacked FLGs 
can be decomposed into groups of BLG bands, with dif- 
ferent effective masses, plus- for odd layer numbers- a 
pair of SLG bands [11]. Fig. 4(c) plots, as an example, the 
schematic band structure of 3LG in a range of the order 
of the phonon energy of the C mode, E(C), and iden- 
tifies electronic transitions that can couple with the C 
mode. Because the density of states with energy higher 
than E(C) is much larger than that with energy smaller 
than E(C), q is not expected to change significantly from 
BLG to bulk graphite, in agreement with our findings. 
If the Fermi energy, Ei?, is larger than E(C)/2, the reso- 
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nance with the C mode would become weaker, eventually 
leading to the disappearance of the BWF profile. 

As discussed above, FWHM(G) is much larger than 
FWHM(C), due to the much larger EPC and phonon- 
phonon contributions. The EPC dominates FWHM(G), 
and the G peak is always lorentzian. 

In summary, we uncovered the Raman signature of the 
interlayer shear mode of FLG. Graphite is not the only 
layered material. Transition metal dichalcogenides, tran- 
sition metal oxides, and other compounds such as BN, 
Bi2Te3, and Bi2Se3 can also be exfoliated to produce a 
whole range of two-dimensional crystals, that are just be- 
ginning to be investigated [35]. Similar shear modes are 
expected in all these materials, and their detection will 
provide a direct probe of interlayer interactions. 
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